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The thermal evolution of neuron stars depends on the elementary excitations affecting the stellar 
matter. In particular, the low-energy excitations, whose energy is proportional to the transfered 
momentum, can play a major role in the emission and propagation of neutrinos. In this paper, 
we focus on the density modes associated with the proton component in the homogeneous matter 
of the outer core of neutron stars (at density between one and three times the nuclear saturation 
density, where the baryonic constituants are expected to be neutrons and protons). In this region, 
f*"^ , it is predicted that the protons are superconductor. We study the respective roles of the proton 

£N| ■ pairing and Coulomb interaction in determining the properties of the modes associated with the 

proton component. This study is performed in the framework of the Random Phase Approximation, 
generalized in order to describe the response of a superfluid system. The formalism we use ensures 
that the Generalized Ward's Identities are satisfied. An important conclusion of this work is the 
presence of a pseudo-Goldstone mode associated with the proton superconductor in neutron-star 
matter. Indeed, the Goldstone mode, which characterizes a pure superfluid, is suppressed in usual 
superconductors due to the long-range Coulomb interaction, which only allows a plasmon mode. 
However, for the proton component of stellar matter, the Coulomb field is screened by the electrons 
and a pseudo-Goldstone mode occurs, with a velocity increased by the Coulomb interaction. 
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I. INTRODUCTION 



Homogeneous matter below neutron star (NS) crust is expected to have a proton superfluid component. The 
elementary excitations of the matter affect the whole thermodynamics and longterm evolution of the star. Since 
the main components of the matter are neutron, protons, electrons and muons [l|, the spectral properties of these 
excitations can have a complex structure. Collective modes in asymmetric nuclear matter have been studied previously, 
e.g. in Refs. [2 4]. In the astrophysical context, a study of the collective excitations in normal neutron star matter 
on the basis of the relativistic mean field method has been presented in Ref. Ja]. The spectral functions of the 
different components in normal neutron star matter have been calculated in Ref. [a0] on the basis of non-relativistic 
, Random Phase Approximation (RPA) for the nucleonic components and relativistic RPA for the leptonic components. 
Different models for the nuclear effective interaction were considered and a detailed comparison was done between 
some Skyrme forces and a microscopically derived interaction. In this work, we will extend this study to the case of 
superfluid matter. The elementary excitations in superfluid neutron star matter have been studied by several authors, 
often with controversial results 043- One important issue is weither the proton superfluid presents a Goldstone 
■ mode at low momentum, since this can have a strong influence e.g. on neutrino emission or mean free path. To 

address this question, we formulate the general theoretical scheme within conserving approximations 15. Il6|. which 
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d . guarantee current conservation and the fulfilment of the related Generalized Ward's Identities (GWI) 

As it is well known, a neutral superfluid must present a Goldstone mode at low momentum. For a single charged 
superfluid, the Coulomb interaction suppresses the mode, which is replaced by the plasmon mode. In neutron star 
matter, the physical situation is complicated by the multi-component structure. The plasmon mode is in fact mainly 
an electron excitation; furthermore the nuclear interaction couples neutron and proton excitation modes. We will 
focus the study only on the proton and electron components, leaving the complete treatment to future works. The 
reason of this choice is that the role of Coulomb interaction is crucial in this case, while the coupling with neutrons 
is expected to be weak [|| . What is left out from the treatment is the possible role of the entrainment [HI, when 
both neutrons and protons are simultaneously superfluid. However, below the crust, neutron superfluidity rapidly 
disappears at increasing density [27l |: therefore we will assume that only protons are superfluid. The aim of the paper 
is to have a qualitative picture of the excitation spectrum and strength function associated with the proton superfluid. 
For simplicity, in the numerical applications, we will consider the effective mass equal to the bare one, and the nuclear 
proton-proton interaction will be treated only schematically. 

The plan of the paper is the following. Section [H] is devoted to the formalism within the conserving approximation 
scheme: the basic equations are derived and their properties discussed. In Section HITl we consider the case of a pure 
pairing interaction. This case is well discussed in the literature, but it is useful to present it before the subsequent 
analysis. In a nuclear system, the situation corresponds to strong pairing, where the pairing gap can be a substantial 



2 



fraction of the Fermi energy. In Section IIV1 we consider the case where both pairing and Coulomb interaction are 
present, and we discuss the interplay of these two interactions in determining the overall structure of the excitation 
spectrum. In Section [V] we present the numerical results for the spectral functions in different physical situations. 
Conclusions are drawn in Section [VIII 



II. FORMALISM 



Let us first consider a non-supcrfluid system. We will concentrate on the density fluctuations. The correlation 
function A is defined by 

A(l,2;3) =<T(p(3)^(l)V(2)>, (1) 

where i = (rj, ti, <7j, Tj) stands for coordinate, time, spin and isospin variables for single particle states; p(3) = 
■0^(3)^(3) is the spin-isospin density operator. According to the conserving approximation scheme [l6j . A satisfies 
the integral equation: 

A(12;3) = A (12;3) + A (12;T 7 2 7 )V(T 7 2 7 ;45)A(45; 3) (2) 

where Ao is the free correlation function and V is the effective particle-hole interaction. A bar over a symbol i indicates 
integration and summation over the corresponding set of variables. Ao is defined by: 

A (12;3) = -G(13)G(32). (3) 

i 

The quantity G is the single particle Green's function: 

G(12) = -i< T{1>(l)tf(2)} >, (4) 

where are the creation and annihilation operators for the considered particles. The effective particle-hole 

interaction V is the key quantity. If the approximation has to be conserving, V must be expressed by the functional 
derivative: 

V(l*45> . (5) 

where £ is the single particle self-energy, defined according the Dyson's Equation 

G _1 (12) = Go 1 (12)-J7(12)-S(12). (6) 

Go is the Green's function for non-interacting particles and U is a possible external single particle potential, which 
can be local in space and time. In Eq. ([5]), the functional derivative is performed considering the Green's function 
as a functional of the external potential. If the derivative is taken at U = 0, then A describes the linear response 
of the system. In this framework, if we express the self-energy in terms of the Green's function according to some 
approximate scheme, the coupled equations @ and © for A and G (or S) define the corresponding conserving 
approximations. It has to be stressed that G and S must be calculated self-consistently. 

The simplest approximations for £ arc the Hartree or Hartree-Fock ones, where the self-energy is a linear functional 
of the Green's function. In this case, the particle-hole interaction V in Eq. (JTJ) is just the bare nucleon-nucleon (NN) 
interaction. However, in general, the bare nuclear interaction cannot be used to calculate the linear response, because 
of the hard core typical of the NN interaction. One can then introduce effective interactions, derived for instance 
from Skyrme forces or from a microscopic procedure, assuming that the interaction can be considered unaffected by 
the nuclear dynamics at the linear response level [6[. A particular case is the Coulomb interaction, which can be 
introduced without modifications, since it contains only a soft core. 

The excitations of the proton-electron system in neutron star conditions were extensively studied in Refs. [H, 0], 
in the absence of pairing effects. It was shown that the electron screening suppresses the proton plasma excitation, 
which is then replaced by a sound-like mode, where the proton and electron components move in phase and are both 
strongly excited. On the contrary, the electron plasma excitation is almost unaffected and remains nearly a pure 
electron mode. 

The inclusion of superfluidity can be formally achieved by including in the above scheme a further discrete variable a, 
labeling the destruction (a = I) and creation (a = — I) operators. The collective index becomes i = (r i7 ti, ai, Tj, Q!j) = 
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(xi,oii). Some care must be used in generalizing several relationships from the normal to the supcrfluid case. In 
particular, the equation defining the inverse Green's function G(l,2) _1 must be written: 



f G- 1 (l,x 3 a 3 )G(x 3 -a 3 ,2) = 6(1-2) 

J rsoraafa 



(7) 



The effective interaction now includes both a particle-hole component V p h and a particle-particle pairing component 
£/p a i r . As a result, the single particle self-energy £(1,2) contains a normal part, and an anomalous part, E a , 
corresponding respectively to a\ ^ ai and a\ = In the Hartree-Fock (mean field) approximation, they can be 
written: 



£"(1,2) 
£ a (l,2) 



JcckJ^i - t 2 )5(r 1 - r 2 )S a 



u ph (ri - ri)G(x'+ - 1, x[ 1)5^ - h) 



-iS(h - t2)S ai a 2 



< XlXilUpairlx'xx'z >A G^tfVjai, r 2 t' 1 <j' 2 a2)5(t' 1 - tl) 



(8) 
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The superscript + indicates that the time variable has been shifted by a positive infinitesimal amount to ensure the 
correct time ordering. The index A indicates the anti-symmetrization of the interaction matrix element. Having in 
mind the Coulomb potential case, we have assumed that the particle-hole component v p h is a density-density local 
interaction, and we have neglected the exchange term; the generalization to a non-local interaction is straightforward. 
The interaction V to be used in the RPA equation can be obtained by functional differentiation, according to Eq. (|5|). 
If, in performing the functional derivative of S, we assume that v p h and t/ pa ir do not depend on G, the approximation 
is still conserving and the procedure is straightforward. More details are given in Appendix [Aj In the following, we 
will consider that pairing occurs in the 1 Sq channel, so that the corresponding effective interaction will be coupled to 
zero total spin. 

Let us now use the momentum representation, obtained by a usual Fourier transform. For simplicity, we assume 
the pairing interaction to have a BCS structure, with a constant pairing gap A and a cutoff (which corresponds to a 
contact interaction in coordinate representation) . Then the generalized single particle Green's functions in momentum 
space are [20| : 

G 12 (q) = ai[5 ai - a3 6 aittr3 G(q)+cTi5 

01,02 $<7\ ,— CT2 

where q is the energy- momentum four- vector (q, co) , and: 
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where q = |q| and m is the particle mass. The equations become algebraic equations directly for the generalized 
polarization function II(q, uj): 



exp[i(qr — qV)]A(r, r; r') = S(q — g')II(q, uj) 



(10) 



where we have used the four- vector notation qr = q • r — cut. The spin and alpha indices are implicit. 

Following the procedure outlined above, we obtain a system of four coupled equations, of which only three are 
independent (one being a linear combination of the others). The generalized RPA equations for the proton system 
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can then be written as: 
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where the X quantities correspond to different components of the Lindhard function, when generalized to the superfiuid 
case: they will be presented in the following. J7 pa i r , the strength of the pairing interaction, is defined to be positive, 
and it is linked to the pairing gap A by the gap equation: 



d 3 k A 

(2tt) 3 2£ k 



d 3 k 

(M 



To express the generalized Lindhard functions, we have introduced the notation: 



l -[X^ G {q)+X^ G {- q )]±X FF {l) 



XT = X G h G (q)±X FF (Q) 
X^ F = X GF (q) ± X GF (-q) 



XGGC = l[X G P G (-q) 



yPP 
GG 



(<?)] 



where the different terms are the following four-dimensional integrals 

dk 



X G h G {q) 

Xl P G {- q ) 
X GF (q) 

X GF (-q) 
X FF (q) 



dk 

— G(k)G(-k + q) 

Jt_ G (k)G(-k-q) 
dk 

^G(k)F(k + q ) 



G(k)G(k + q) ; X G h G {-q)=X G %{q) 



rph 



X FF (-q) = X FF (q) 



The explicit expressions of X^P, X^ and X GF are given in Appendix [C] 

If the alpha indices are explicitly indicated, the polarization tensor can be written: 

n(q,w) = n aia < ;a2C> /(q,a;) 

and the three components of the polarization tensor appearing in Eq. arc defined as 

1 
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(23) 

(24) 
(25) 



where the values of the variables a and /3 are generic (they are chosen according to the polarization tensor components 
that have to be calculated) . The index S in Eq. ([TTj) specifies that scalar (zero total spin) excitations are considered. 
The BCS polarization tensor appears on the right hand side of Eq. (fTTj). It can be expressed in terms of the same 
quantities X given by Eqs. (|13H15|) (sec Appendix [UJ. 

In the so-called constant level density approximation, i.e. exact particle-hole symmetry, the off-diagonal matrix 
elements of Eq. (JTTJ) , which couple the first equation with the other two, are of order q 2 , and they contribute to the 
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determinant only through terms of order q 4 . As a result, to order q 2 , the first equation is decoupled from the other 
two, and actually can be neglected. 

Finally one must take into account the Coulomb coupling between protons and electrons. Including proton and 
electron components, the RPA polarization tensor forms a system of three coupled equations: 
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(26) 



where X e is the relativistic Lindhard function for electrons (2l|, calculated in the Vlasov limit, and n( ee ) is the 
corresponding part of the polarization tensor involving the electron (density) component. These are the equations to 
be solved to get the proton and electron strength functions in the presence of pairing correlations. 

To extract the proton density-density polarization tensor n_ii ; _n, it can be of interest to express the electron 
polarization tensor in terms of the other components, and to obtain a reduced set of equations for Hg + ^ and Hg\ 
One gets: 
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where the screened Coulomb interaction is given by 



v c /(l - 2X e v c ) 



(27) 



(28) 



This effective interaction is energy dependent. If one approximates the polarization function X e by taking the zero 

~(s) 

frequency limit at fixed momentum, v c becomes a static screened Coulomb interaction Vc between protons: 



(q) = 47re 2 /(q 2 + q 2 ) 



(29) 



The screening wave number is given by q 2 = 3(Wp/wFe) 2 , where is the electron plasma frequency and vp e is the 
corresponding Fermi velocity. We will refer to this approximation as the " static approximation" . 



III. PAIRING CORRELATION ONLY 



The simplest situation is when only the pairing interaction is included. For the considered case of Eq. (|26[) this is 
formally equivalent to put v c = 0, and it corresponds e.g. to an uncharged superfluid. The corresponding expression 

for Tl^g^ can be easily obtained by solving the remaining two by two algebraic system. Replacing the r.h.s. of Eq. (|26[) 
by its expression in terms of the quantities X, given in Appendix [C] we obtain: 



IT 
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ph 



/(i-xf[/ pair ) 



(30) 



One can verify, as shown in Appendix [Cl that the numerator of this expression vanishes at q = for any non-zero value 
of uj. This implies that the response function and the corresponding strength function, which are even functions of q, 
are proportional to q 2 for small q for any non-zero value of oj. This is a consequence of the conserving approximation 
we are following, which guarantees the conservation of current: the continuity equation implies this property for the 
density-density response function (see Appendix This result is in agreement with Ref. (22j . where the particular 
expression for the vertex function from Ref. [Toj was adopted; however this expression was valid only for oj <C 2A. 
In any case this RPA property of current conservation is a well known result in the theory of superconductors [l7| . 
where only pairing interaction is present, and it is also a consequence of gauge invariance. 

According to a general theorem on symmetry breaking, the strength function for u < 2A is characterized by the 
presence of a "Goldstone mode", i.e. a phonon-like excitation with an energy proportional to q for small enough q. 
It is well known [l7j that in this limit, and in the weak coupling approximation, the velocity of the Goldstone mode 
is vf/V3, where vp is the Fermi velocity. We see on Fig. [T] that similar features are obtained in our case. This figure 
represents the excitation spectrum of a proton superfluid, in the absence of Coulomb interaction, for a proton density 
corresponding to neutron star matter at total baryonic density equal to the nuclear saturation density po = 0.16 
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FIG. 1: Proton excitation spectrum with pairing interaction only. The proton density corresponds to the case of neutron-star 
matter at nuclear saturation density, with a proton fraction fixed by /3 equilibrium, which means a proton Fermi momentum 
kF P = 0.56 fm _1 . The chosen pairing gap is A = 0.5 MeV. 
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TABLE I: Density dependence of the proton fraction, Fermi momenta and electron plasmon frequencies in neutron-star matter 
conditions. The values given for the proton fraction Y p are obtained by microscopic calulation. 

fm -3 . The proton fraction (determined by f} equilibrium in stellar matter) is taken from Ref. 0: it is reported for 
convenience in Table |U Here we assume for illustration a typical value A = 0.5 MeV. The exact value of A is not 
known, since the microscopic many-body theory of pairing in nuclear matter is not a settled issue (24|. The energy 
of the excitation branches are defined as the zeroes of the determinant of the real part of the matrix on the left hand 
side of Eq. (|26|) , for each given value of the momentum q. 

For small q one observes a linear behavior of the Goldstone mode energy; a numerical measurement of the slope 
shows that the corresponding velocity is indeed very close to vf / \/3. This velocity can also be calculated according 
the equation: 

1 - C/pah-^f = (31) 

which corresponds to the vanishing of the determinant of the matrix on the left hand side of Eq. (|26|). Using an 
expansion of the quantity X^ up to order q 2 , and taking the limit A <C E-p , we recover from Eq. pip the standard 
Goldstone velocity vp/y/3. Let us note that the analytical result depends only weakly on the value of A/Ep: this 
explains why, in our strong coupling case, we obtain a Goldstone velocity similar to that of the weak coupling limit. 
More details on the analytic expression of the Goldstone mode velocity are given in Appendix ID] 

According to Eq. ([3171) . the strength function has a delta singularity at the energy of the Goldstone mode, since it 
is undamped in the considered limit of pairing interaction only. One can see on Fig. [T] that the linear trend is not any 
more valid when the energy is approaching 2A. Actually the mode can exist only below 2A, since otherwise it could 
decay in two quasi-particle states and it would be completely damped. Above 2A one observes another branch, that 
starts at 2A for q = 0. This excitation can be identified with a "pair breaking" mode. A low-momentum development 
of Eq. (j3"Tj) also gives an analytic expression for the energy of this mode, and shows that it is proportional to q 2 . 
However, this quadratic dependence disappears for larger values of q. Both branches tend to an energy value close to 
2A: this could be expected, since the quantity X 1 ^ has a logarithmic singularity at uj =2A for all reasonable values 
of the momentum q, typically for values smaller than the Fermi momentum. Thus, the two branches nearly touch at 
a given momentum, where they both disappear. 
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FIG. 2: Proton excitation spectrum with pairing and Coulomb interactions, ignoring the electron component. The proton 
density corresponds to the case of neutron-star matter at nuclear saturation density, with a proton fraction fixed by /3 equi- 
librium, which means a proton Fermi momentum ltFp = 0.56 fm _ . The pairing gap takes different values, indicated by the 
label: A = 0, 0.25, 0.5 and 1 MeV. 

IV. DEATH AND RESURRECTION OF THE (PSEUDO-)GOLDSTONE MODE. 

It is a classical result [23| that in a charged superconductor the Goldstone mode cannot be present: the standard 
theorem on broken symmetry (the gauge invariance in the present case) is not any more valid if the particles interact 
also by a long range interaction like the Coulomb one. The mode is then replaced by a plasmon mode, whose frequency 
is only slightly modified with respect to the case of a normal charged fluid. The plasmon branch starts from a non zero 
energy at q = (the plasmon energy) and varies slowly with the momentum. This is what we obtain if we consider 
the first two lines of Eq. (|26l) . including the Coulomb interaction v c but ignoring the electron component: then we 
deal with a charged superconductor, the proton liquid. In this case, the determinant of the matrix on the left hand 
side of Eq. (|26[) vanishes at non-zero energy for q = 0, since the factor in front of v c actually vanishes, as it happens 
in a normal charged fluid [f|. 

The full structure of the excitation spectrum is illustrated in Fig. [5J for the same proton density as in the previous 
section, and for different possible values of the pairing gap. When the gap A is much smaller than the plasmon energy, 
the structure of the branches resembles the "thumb like" shape typical of the charged normal Fermi liquid (3. l2lll25ll26| . 
In this case, the lower branch is over-damped: it does not correspond to an actual excitation. The upper branch 
corresponds to the plasmon mode, which is mainly undamped up to a momentum nearby the end of the "thumb", 
above which no excitation mode exists. In the present superfluid case, the lower branch must start above the forbidden 
energy region uj < 2A. At increasing value of A, the spectrum changes considerably, although the position of the 
genuine plasmon excitation remains nearly unchanged: for a large enough value of the gap, the momentum dependence 
of the plasmon energy deviates more and more from the quadratic form typical of a charged Fermi liquid, and the 
whole spectrum shrinks. When 2A reaches a value close to the plasmon energy uj p of the normal proton fluid, no 
excitation is possible any more. 

On the other hand, in a physical situation such as neutron-star matter, we have to take into account the presence 
of the electrons, which has a crucial impact on the excitation modes of the medium. The structure of the spectrum 
changes radically when the electron component is introduced. Indeed, electrons are much faster than protons, and are 
able to screen the proton-proton Coulomb interaction at all considered frequencies. As a result, the effective proton- 
proton interaction is of finite range, and a sound-like branch appears again. This effect was extensively discussed 
in [f| 0] in the case of non-superfluid matter. In the present situation, below 2A, the sound-like branch associated 
with the proton component can be considered as a pseudo- Goldstone mode, which is determined both by pairing and 
screened Coulomb interaction. 

The structure of the excitation spectrum is illustrated in Fig. [31 for different values of the pairing gap. In this 
figure, we report the branches where the strength of the proton component is large. The plasmon mode, which is 
mainly electronic, does not appear: it occurs at higher energy. The first panel shows the two branches associated 
with a normal proton fluid, with no pairing. The upper branch corresponds to a sound mode, which presents some 
damping, and whose velocity is determined by the screened Coulomb interaction. The lower branch is over-damped 
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FIG. 3: Proton excitation spectrum with pairing and Coulomb interactions, including the electron component. The proton 
density corresponds to the case of neutron-star matter at nuclear saturation density, with a proton fraction fixed by /3 equi- 
librium, which means a proton Fermi momentum kp p = 0.56 fm _ . The pairing gap takes different values, indicated by the 
label: A = 0, 0.25, 0.5 and f MeV. 



and does not correspond to a true excitation mode. The two branches join and stop at a cutoff momentum. When 
the proton pairing is introduced, the spectrum is still composed of two branches, but a clear distinction appears 
between two momentum domains. For low enough values of q, one single branch is present below 2A, with an energy 
approximately proportional to the momentum: this is the pseudo-Goldstone mode. In this domain, the upper branch 
is located above 2A: it corresponds to the pair-breaking mode. This region of the spectrum is quite similar to the case 
of pure pairing presented in the previous section, although the pseudo-Goldstone velocity is different from the pure 
Goldstone one, as will be discussed in the following. For increasing values of q, these two branches get closer, until a 
given momentum where they undergo a quasi-crossing. Afterwards, the upper branch becomes similar to the sound 
branch of the normal fluid, and like for the normal fluid, the lower branch is over-damped and does not correspond 
to a true excitation. One can see that the part of the spectrum above 2A tends to shrink as the value of the gap 
increases. For a large enough gap, it actually disappears and only the pseudo-Goldstone and the pair-breaking modes 
are present, as illustrated in the last panel of Fig. [31 at A = 1 MeV. 

As noticed in the previous section, if only the pairing interaction is considered, a Goldstone mode is present below 
2A, with a velocity vq that in the weak coupling limit is equal to vpp/^/3. The pseudo-Goldstone mode that is 
obtained in the presence of a screened Coulomb interaction has a much higher velocity, about three times larger. 

The excitation spectrum in the region of small values of the momentum q can be studied analytically. As it is 
shown in Appendix [Dl the expansions of the matrix elements appearing in Eq. (|26[) give an explicit formula for the 
velocity vpq of the pseudo-Goldstone mode. If the small cut-off dependent terms arc neglected, in the weak coupling 
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FIG. 4: Velocity of the Goldstone mode as a function of proton density, vg : velocity of the Goldstone mode in the pure pairing 
case (purepair). vpg '■ velocity of the pseudo-Goldstone in the presence of proton pairing and screened Coulomb interaction, 
calculated in the static approximation (static) or with an exact treatment (exact). 



limit and static approximation for the electron polarization function one gets a quite simple formula: 



y PG 



Np 



(A/E Fp 



(32) 



where N p and N e are the proton and electron level densities. Since the protons can be considered non-rclativistic and 
the electrons ultra-relativistic, one easily obtains 



Np 
N P 



m p c 



Fp 



(33) 



where m p is the proton mass. At the considered density this ratio is about 8, which gives the factor 3 between 



the pseudo-Goldstone and Goldstone mode velocities. Equations (|32|) and (|33|) show that the ratio between the two 
velocities decreases with increasing density. 

The static approximation for electron polarization function is expected to be quite accurate at low momentum and 
energy. Figure Q illustrates the relevance of this approximation for the calculation of Vpq. The numerical calculations 
are in close agreement with the analytical results and the agreement is present at all densities. Furthermore, we have 
verified that vpq/vq has a rather weak dependence on the value of the gap. 



V. THE STRENGTH FUNCTIONS. 



The density strength function 5„(q, w) is the key quantity for many physical properties which characterize the 
behavior and dynamics of neutron star cores, like neutrino scattering, neutrino emissions, specific heat and cooling. 
In the present notation, it is given by: 

S„(q >W ) = (34) 

where 3 indicates the imaginary part. The function 1$?^ is also called the vector response function, in contrast with 
the axial vector response function which corresponds to the spin-density strength function. The latter, which is also a 
key quantity, can be studied along the same lines; however, we will concentrate on the vector response function. The 
strength function in the absence of pairing interaction was studied in Ref. @ . The pairing correlations have different 
effects, that we analyze now in detail. 

The behavior of the strength function reflects the structure of the excitation spectra obtained in the previous 
sections. The different peaks that appear correspond to the branches in the spectrum. However the position of the 
peaks in the strength functions can be slightly shifted with respect to the corresponding branches. This is due to the 
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FIG. 5: Proton spectral function when only the pairing interaction is included. The proton Fermi momentum is kF P = 0.56 
fm" 1 . 
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FIG. 6: Proton spectral function when pairing and Coulomb interactions are included, without electron screening. The labels 
indicate the different values of the pairing gap in MeV. The momentum is q = 0.025 fm -1 . The proton Fermi momentum is 
k Fp = 0.56 fm" 1 . 

imaginary part of the response function, i.e. to the damping of the excitations, which moves the position of the mode 
with respect to the branch extracted from the real part only of the RPA matrix. Furthermore, it can happen that a 
branch is over-damped and does not correspond to a real mode, like it is the case for the branch below the normal 
sound mode in Fig. |3l 

For the sake of illustration and analysis, let us first consider, as in Section IIII1 the case of a proton system with 
only the pairing interaction, and no Coulomb effect. At not too high momentum, two branches are present in Fig. [1] 
In the corresponding strength function, these two branches give rise to two structures, as shown in Fig. [5j where the 
proton strength function at q = 0.04 fm -1 is reported. The narrow peak below 2A corresponds to the Goldstone 
mode, which is undamped, i.e. it corresponds to a Dirac-delta singularity (numerically one gets just a sharp peak). 
Exactly at 2A, the spectral function jumps: this corresponds to the opening of the two quasi-particlc channel, i.e. 
to the possibility of breaking a Cooper pair. Above 2A, we can observe a structure that corresponds to the pair 
breanking. It is not sharp because it is immersed in the continuum. 

We now introduce the Coulomb interaction. Like in Section HV1 we first consider the electrons as a static homoge- 
neous background, which has no role in the excitation spectrum of the protons. In this case, the strength function 
is dominated by the proton plasma excitation above the forbidden region. According to Fig. [3J the proton plasmon 
mode (upper branch) remains essentially unchanged at low momentum, but the spectrum shrinks at increasing A and 
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eventually disappears. The corresponding spectral function is shown in Fig. [6] for q = 0.025 fm -3 , and for different 
values of the gap. The plasmon excitation is well pronounced when the pairing gap is small, since then we are close 
to the normal fluid case. As A increases, this mode is more and more damped and shifted. We can notice that 
the damping is strong even when the position of the excitation along the branch is well below the endpoint of the 
spectrum (the tip of the " thumb" ) . As regards the lower branch, it is always over-damped, as in the normal system: 
no strength is associated with it. We will not develop further the analysis of this case, since the role of electrons has 
to be introduced in order to discuss the realistic physical situation in neutron stars. 

Let us now consider the strength function when the electron component is treated dynamically, as in a realistic 
situation. This case is illustrated in Fig. [Jj for q = 0.025 fm -1 and for different values of the pairing gap. When 
the pairing gap is very small (upper left panel), the system is close to the normal one: the spectral function is 
dominated by the sound mode, well above the so called forbidden region u < 2 A. The mode appears as a relatively 
sharp peak, where the electron strength is comparable with the proton strength. This indicates the screening effect 
of the electrons, which are much faster than protons and can easily follow their motion. The electron screening is 
responsible for the transformation of the proton plasmon mode into the sound mode when the electron component is 
treated dynamically. The position of the peak can be inferred from the upper branch shown in the first panel of Fig. [3] 
The width of the sound mode is determined by the proton-electron coupling, since its position is above the proton 
particle-hole continuum, so that direct Landau damping can not occur. As usual, the lower branch is over-damped 
and no strength is associated with it. 

In the upper right panel of Fig. [JJ the pairing gap A is increased to 0.2 MeV. Then, the sound mode falls closer to 
the forbidden region, but still above it, as can be infered from the upper right panel of Fig. [3] The mode is strongly 
damped by the coupling with the two quasi-particle continuum: the Landau damping is now active. 

In the lower left panel of Fig. A is further increased to 0.5 MeV. A sharp peak appears again at low energy, but 
now it is located below 2A: this is the pscudo-Goldstone mode. In contrast with the Goldstone mode of the pure 
pairing case, wich corresponds to a Dirac singularity, the pscudo-Goldstone mode is slightly damped: this damping 
is due to the coupling with electrons, since no proton pair breaking can occur in this region. Above 2A, a broad 
structure is present, corresponding to the pair-breaking mode. Although its strength is spread over a large energy 
interval, its total strength turns out comparable with the one of the pseudo- Goldstone mode. Further increase of the 
pairing gap (lower right panel) confirms this trend of the strength function. 

In the last panel, one observes also an abrupt decrease of the strength function at a given energy (about 5 MeV): 
this is due to the behavior of the electron free polarization function X e (relativistic Lindhard function) involved in 
Eq. (1261) . Indeed, for q <C kp p , the imaginary part of X e is proportional to q, and contributes to the imaginary part 
of the particle- hole polarization propagator W h , In the ultra-relativistic limit [2l|, this contribution suddenly goes 
to zero at E = ficq: this causes the sharp drop of the strength function. A similar effect, but less pronounced, is 
observed in the lower left panel. This also shows the mutual influence between electron and proton components in 
the excitation modes and the corresponding spectral functions. 

The electron strength function, besides following closely the proton strength function at lower energies, displays 
a sharp plasmon mode at very low values of the pairing gap. However, as soon as the gap increases, the electron 
plasmon mode is strongly damped. Since the proton Fermi function is smoothed by the pairing interaction, we 
interpret this effect as due to the direct coupling of the electrons with the two proton quasi-particle excitations, whose 
continuum can extend up to the position of the electron plasmon mode. Notice that in this region the position of the 
electron plasmon mode is well above the electron particle-hole continuum and no direct Landau damping is possible. 
The electron plasmon mode actually disappears at values of the gap of few hundreds of keV. Such a strong effect is 
probably overestimated in our calculations, where a constant gap has been assumed up to twice the Fermi energy. If 
the pairing gap was assumed to be strongly concentrated around the Fermi momentum, the coupling and damping 
effects would be reduced. However, the results clearly display the qualitative behavior of the electron strength function 
as the pairing gap is increased. 

In Fig. [8] is illustrated the development of the proton spectral function, which is mainly concentrated at low energy, 
as the momentum increases, for a fixed value of the gap A = 0.5 MeV. At low momentum (upper left panel), the main 
mode of excitation is the pseudo- Goldstone mode below 2A, together with a well defined structure corresponding 
to the pair- breaking mode, as discussed above. As the momentum increases (upper right panel), the main peak 
moves close to 2A, but still keeping below, while the pair-breaking mode is much more spread. Notice that still the 
electron strength function follows closely the proton one. Increasing further the momentum (lower left panel), the 
pseudo- Goldstone mode turns into the very damped sound mode above 2A, with a residual strength at 2A. Finally, 
at higher momentum (lower right panel) , the sound mode develops further and becomes slightly more pronounced. As 
the momentum increases, the electron strength is reduced with respect to the proton one, since at higher frequency 
the electron screening becomes less effective. 

Another feature of these spectral functions is the presence of two kinks. The lowest in energy is trivial: it corresponds 
to the opening at 2A of the two quasi-particle channels, i.e. to the breaking of the Cooper pairs. The second one, 
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FIG. 7: Proton and electron spectral functions, with proton pairing and Coulomb interaction, for a fixed value of the momentum 
q = 0.025 fm _1 and different values of the pairing gap A. The proton Fermi momentum is kF P = 0.56 fm~ . 



indicated by an arrow in the lower panels of Fig. is less obvious. It can be understood from the kink that is present 
in the two quasi-particlc density of states. One finds that it occurs at the energy ( q <C kp p ): 



kink 



2\ e n E 



A 2 



(35) 



With a momentum dependent pairing gap, the kink would be smoothed, but in any case a sudden change of slope is 
expected to occur. 



VI. INTRODUCING THE NUCLEAR INTERACTION 



The effective nuclear particle-hole interaction in the very asymmetric matter present in neutron stars is not well 
known. At small momenta, as in all Fermi liquid, it can be characterized by a set of Landau parameters. In 
the language of Landau theory, the density-density interaction F(q) of two quasi- particles at the Fermi surface is 
the particle-hole interaction to be used in the RPA equations at low momenta. It can be expanded in Legendre 
polynomials as a function of the "Landau angle" between the momenta of the two quasi-particles: 

F(q) = Fflki - k a |) = ^F L P L (cos0) (36) 

L 

where |kx| = |k 2 | = kp. For asymmetric matter, one has three interaction terms, corresponding to neutron-neutron, 
neutron- proton and proton-proton channels. The term Fq, corresponding to L = 0, is usually dominant. In Ref. 0, 
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FIG. 8: Proton and electron spectral functions, with proton pairing and Coulomb interaction, for a fixed value of the gap 
A = 0.5 MeV and different values of the transfered momentum q. The proton Fermi momentum is kF P = 0.56 fm -1 . The 
arrows point to the kink discussed in the text in relation with Eq. (|35[) . 



it was shown that it can be obtained from the single particle potential at the Fermi surface, as obtained from the 
Brueckner theory. If we restrict to the proton-proton channel, the interaction matrix element V^ c at saturation total 
baryon density turns out @ to be repulsive and about 311.5 MeV fm 3 . It corresponds to a Landau parameter 

F « = V ™^~ 0.426. (37) 

The relevance of such an interaction can be estimated by comparing it with the screened Coulomb interaction of 
Eq. (j2"9")l . Even at q = 0.1 fin -1 , the nuclear interaction is less than one fourth of the Coulomb one: therefore, only 
minor corrections are expected to arise from the nuclear proton-proton interaction. In fact, one finds that only near 
the endpoints of the spectrum of Fig. [3] some corrections to the energies occur, not larger than 10%. As a further 
illustration, we show in Fig. [9] the proton and electron spectral functions with and without the nuclear proton-proton 
interaction. One can see that the nuclear interaction, being repulsive, slightly shifts the proton spectral function to 
higher energies. In any case the effect is only marginal. 

The effective nuclear interaction can be also calculated from Skyrme forces. In a previous paper we have 
compared the effective forces calculated from the microscopic approach with the ones coming from a set of Skyrme 
forces. Despite some qualitative similarities, the effective forces are quantitatively quite different. The calculation of 
the effective forces for such highly asymmetric matter is a challenge both for the microscopic and the phenomenological 
approaches. 

As for the effective mass, the same microscopic Brueckner calculations gives a proton effective mass to* at the Fermi 
surface of about 0.89 to, being to the bare mass. It is well known that the inclusion of dispersive dynamical effects 
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FIG. 9: Proton and electron spectral functions, with proton pairing and Coulomb interaction (thin lines); the thick lines 
represent these spectral function when the nuclear proton-proton interaction is included. The pairing gap is A = 0.5 MeV and 
the proton Fermi momentum is k.F P — 0.56 fm _1 . 



enhances slightly the effective mass, that is then expected to be close to m. In conclusion, also the introduction of a 
realistic proton effective mass can only marginally affect the results presented in this paper. 



VII. SUMMARY AND PROSPECTS 



We have presented a study of the elementary excitations in the homogeneous asymmetric nuclear matter that is 
present in neutron stars, focusing the analysis on the proton component. The protons have been assumed to be 
superfluid, and we have considered the pairing gap A as a parameter taking values up to 1 MeV. We have shown how 
the Goldstone mode that is present in a pure uncharged superfluid below 2A becomes a pseudo-Goldstonc mode due 
to the Coulomb interaction. This mode has a velocity about three times larger than the original Goldstone mode. We 
have calculated the spectral functions at different momentum: at increasing momentum, the pseudo-Goldstonc mode 
merges, through a quasi-crossing mechanism, into the sound mode that characterizes the normal system. Besides these 
modes, the pair-breaking mode is clearly apparent in the spectral functions. However it has a quite broad structure 
with a modest strength. Finally, we have considered the influence of the nuclear proton-proton interaction. It turns 
out that it has only a weak effects on the overall structure of the spectrum and of the spectral functions. 

The main items that still have to be analyzed are the introduction of the neutron components and the study of the 
axial-vector channel, with the ultimate aim of calculating the neutrino emission rate and mean free path. Work is in 
progress in this direction. 



Appendix A: Conserving approximations 

In this appendix, we discuss the conserving approximation for the response function. This approach was developed 
a time long ago by G. Baym [l6| for a normal system, with a central particle-particle interaction. The theory can 
be readily generalized to a more complex interaction and to superfluid systems. The starting point is the equation 
of motion for the single-particle Green's function G, which involves the two-particle Green's function G2. Then it is 
possible to show that if one approximates Gi by a functional of G. which satisfies a definite functional condition, the 
average particle density < n(r, t) > and current < j(r, t) >, calculated in terms of the single-particle Green's function, 
satisfy the local conservation law: 

^-<n{r,t)> +V<j(r,i)>= (Al) 
at 

This result holds true if (i) the single particle Green's function satisfies self-consistently the equation of motion 
with the approximate G2 and (ii) the interaction is local. The main reason of this result is the fact that, under the 
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functional condition on G2, the contribution of the interaction term to the current vanishes: the conservation law 
expressed by Eq. (|A1|) then follows trivially from the kinetic term contribution. 

For simplicity, we illustrate this point in the case of a pairing interaction in the mean field (BCS) approximation 
for G2 in the equation of motion: this approximation is known to be conserving. Let us take a zero-range pairing 
interaction, which is the one we have actually used in the present work: 

^pair = -\ U pail 5w J d 3 r^(v,a)^(r,-a)^(r,-a')i;(v,a')(-l)^- a -' 7 '> (A2) 

where spin-zero pairing has been assumed and U pa \ r (> 0) is the corresponding strength. The relevant component of 
the Green's function G(l, 1') corresponds to ai = 1 and = — 1. This component will be indicated by G{x\,x' l ). 
The equation of motion for G can be written: 



,-G(x 1 ,x 1 ) + - 

+1 U paiT E a , < T(^(r-a u t)yj(r- a 1 , t)^(va' , t)^ {v' a^t')) > (_i)d— 1-') 



ift^-G(a;i,xl) + ^-V?G(xi,ii) = S(t - t')6(v - r') (A3) 



where x\ = (r,a\,t) and x\ — (r',ai,t'). The analogous equation for the derivative on the time t' reads: 

m^G{x 1 ,x' 1 )-^V 2 r G(x 1 ,x' 1 ) = S(t-t')8(r-r') (A4) 

+^ U pail X„ <T(^(r* 1 ,t)^(r'a,t')^(r' -<r,t')ip(r' -^t')) > (-1)( 1+ ^' 

Summing up the two equations of motion, performing the summation over the spin 01, and finally putting t — t' and 
r = r', the first two terms on the right hand side give the conservation law (|A1|) . Let us remind that the density and 
current can be written in terms of the single particle Green's function: 

<n(r,t)> = < ^ f (rai,t)tp(ra 1 ,t) > = iE (ri G(x,x + ) (A5) 

<j(r,t)> = A(v r -Vr/) <i> ] (nri.t)^^!, t) > \ r ,= r = E CTl (V r - V r >)G(x, x+)\ r , =r 

where x + = (r, cr\,t + e), with e a positive infinitesimal quantity that enables to fix the correct time ordering of the 
operators in the equal-time limit, according to the Green's function definition. 

In the mean-field approximation, the interaction contribution on the right hand side of Eqs. (|A3[) and (|A4[) can be 
factorized. For the first Eq. (|A3I) . one has: 



< T(1»\T,-a U t)1>(T,-</,t)ll)(Ty,t)1>HT f , > 

« <T(V t (r,-a 1 ,t)V t (r' )t 7 1 ,t / )) >< V(r,-</,t)V>(r,a',t) > 
+ <^( r) _ (ri) ^( r) _ cr ' ) t) ><r(^(r,<r / > t)V t (r'.^i.* / )) > 

- < ^(r, -<ri,t)7/;(r,<7', t) >< T(^(r, -a' ,t)^ (r 1 ,a u t')) > (A6) 



A completely analogous factorization applies in the second Eq. (|A4[) . Taking into account the translational and time 
reversal symmetries, which imply: 

<T(^{r-<i 1 ,t)^(r'cr 1 t r ) >= - < T(^(r - a x , t^r'ax, t') > , (A7) 

the two factorized G2 are equal. Since they are multiplied by an opposite phase in the equations of motion, they 
cancel out. The conservation law for the response function follows by using the basic formula for the functional 
derivative fl6j | : 

6 -^fTT =~i< T^(x)^(x')n'( Xl ) > (A8) 
oil (Xi) 

where U{x\) is an external scalar potential coupled linearly to the density n(xi), and n'[x\) = n{x\)— < n(xi) >. 
Indeed, the functional derivation of Eq. (J7J gives 



■1'3 



5G 1 (l,x 3 a 3 ) , r -i,, 8G(x 3 -a 3 ,2) 

G(x 3 - a 3 ,2) + G (l,x 3 <V 



SU v ' ' woo/ S[J 



(A9) 
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By using the chain property of the functional derivative, one gets the integral equation 

A(12;3) = A (12;3) + A (12;T 7 2 7 )V(T^;;45)A(45;3) (A10) 

where a subscript * indicates a sign change of the variable a, e.g. 1* = (ri,<7i, ti, — a\). The function A has the same 
expression as in Eq. Q, provided the Green's functions are generalized to the supcrfluid case. The correlation function 
A(12; 3) and the integral equation (|A10[) are the generalization to the superfluid case of Eqs. ([1} and ([2]), respectively. 
If one puts x\ = X2 and a\ = —a^ — 1, the density-density response function Il( p ' l )(l, 2) can be recovered. In this 
way the conserving property verified by A is directly transferred to the generalized response function, provided the 
Green's function (or self-energy) is calculated self-consistently. In the mean field approximation, which generates 
the RPA approximation for the generalized response function, this means that the self-energy appearing in the RPA 
equation for the response function must be calculated self-consistently within the mean field approximation. For the 
pairing part this means that we have to use the BCS Green's functions. This property of the RPA equation, when 
only pairing interaction appears, was recognized a long time ago [17| . 

If other interactions than pairing arc present, noticeably the particle-hole interaction, the argument can be repeated 
and this RPA property still holds, provided all interactions are local. The mean field corrections to the self-energy in 
the mean field approximation amounts to a shift in the chemical potential, if the exchange term is neglected (this is the 
Hartrce approximation, which is also conserving). If the interaction is non-local, or if we include the exchange term, 
the local continuity equation (|A1[) does not hold any-more. It can be restored in the effective mass approximation, 
where the effective mass is momentum and energy independent. 

It has to be noticed that the conservation law for the vertex function A, which appears in Eq. (IA10[) . is equivalent 
to the generalized Ward's identities [l7j . The latter imply the local conservation of current, which for the polarization 
tensor in momentum-energy representation can be written 

W ngj + Si =1 , 3 qi ngj = o (Aii) 

where the index value i = corresponds to the density component and the values i = 1, 2, 3 correspond to the three 
current components. The (0,0) component II coincides with the polarization function H^p^ of Eq. (|2"6"|) . For an 

isotropic system, the second (current) term is proportional to q 2 for small q. Thus, the conservation law implies that, 
for any value of the energy w, the density-density response function in the limit of small value of q vanishes as q 2 . 



Appendix B: Self-energy and effective interaction 



Let us consider the equation of motion for the single-particle Green's function in the superfluid case. The nucleon- 
nucleon hamiltonian can be expressed in general as: 

H = H + V 
fi 2 V 2 



Ho 



2m 



(Bl) 



where the integrations include summation over spin (for simplicity we do not include the isospin). All time variables 
are equal. The interaction can be an effective one, in which case it is considered static and particle number conserving 
(weak coupling limit). 

The equation of motion for the single-particle Green's function can be derived generalizing the result of Appendix fA"! 



ih^- — = -i 
at\ 



- ai ^^<T{V>(l)V(2)}> 
2m 



<xiX2>\v\x3>xa> > A < T{tl)\ i {x2')iit 1 {x v )ip tl {x3 l )ip{2)} > 

< X V Xy\ v \ x ^ x U >A< T{tpl l (x 1 ,)^l l (x2')lpt 1 ( x i'H{2)} > 



2 

5 ai ,-i 



+ hS(t 1 -t2)S(r 1 -r 2 )6 ai - a2 5 ait(T2 



(B2) 



where t[ = t 2 = £3 = t' 4 = t±, the subscript A indicates anti-symmetrization and the integration includes summation 
over spin. Remind that ^>(2) = tp( x 2) if «2 = 1, and ip(2) = ij)^(x2) if = — 1- In the Hartree or Hartree-Fock 
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approximation, the two-body Green's functions on the right hand side of Eq. (|B2|) can be now factorizcd in the same 
fashion as in Appendix |A"1 This enables to extract the self-energy according to the Dyson's Equation ©: 

£(1,2) = ffife- _L J [2S a2 - ai < xix v \v\x 2 x 2 > >a (G(x2',tf,-l;xi',tl,l)6 ait - 1 -G(xif,tf,-l;x2',tl,l)6 aitl ) 

+S a2 , ai < xix 2 \v\xi>x 2 > >a G(xi/,tf ,ai;£ 2 ',ii ,ai)] ( B3 ) 

If we now assume that the interaction is the sum of a local particle-hole interaction and a local zero-range pairing 
interaction, we recover the expression of Eqs. (|8I9[) for the self-energy. From Eq. (|B3|) . we can derive the effective 
interaction V, to be used in RPA equations: 

J£(l 2) 

V(l, 2; 4, 5) = ^ = ifi(h - h)5a 2 ,-ai5oi4.,-l5a 5l l [< XiX 5 \v\x 2 X4 >A S ai ,-1- < XiX4\v\x 2 X5 >A 5 au i] 

1 

< xix 2 \v\x±X5 >a (B4) 

Again, if we assume that the interaction is the sum of a local particle-hole interaction and a local zero-range pairing 
interaction, the expression simplifies and we can easily derive the RPA equations in the energy-momentum represen- 
tation that we have used in our calculations: see Eqs. (fTTj) and (|26|) . 



Appendix C: Generalized Lindhard function and numerical method 



The relevant functions X appearing in the kernel of the RPA equations, reported in Eqs. (|13H15[) . after the integra- 
tions over energy, can be written: 
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•£|k+ q | + E k -ui - ir) £jk+q| + E k + oj - it] 



Despite the approximation of constant pairing gap, these integrals cannot be calculated analytically and their numer- 
ical evaluation requires particular care. We follow the method of assigning to the quantity r\ a small value, simulating 
an infinitesimal, and then integrating by an adaptive method the remaining two-dimensional integrals. We have used 
the subroutine DT20DQ of the Visual Fortran package. In this way, both real and imaginary parts can be calculated 
in few minutes with a simple PC. We have found that the results are quite stable if we take r\ values between 10~ 3 
and 10 -6 (energies are all calculated in MeV). For production calculations we have used the value rj = 1.5 x 10 ~ 5 . 

Finally, for completeness, we give here the expression of the free polarization functions appearing on the right hand 
side of RPA Eq. (|26[) . In the case we are dealing with, i.e. the density-density response function, which corresponds 
to a = — 1 and /3 = 1 in Eq. (|2l)|) . we have: 
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Using these expressions, for the pure pairing case (i.e. v c — > 0) one gets the expression of Eq. (J3UJ) for the density- 
density component II^'' 1 '. Then, putting q = in the previous expressions for the functions X, one finds: 
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J pair 



to 2 A 2 



d 3 k 



(2tt)3 E k [{2E k f-^ 



= X p _ h (1 - u } 



vPP\ 
pair ^+ ) 



(CI) 
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Together with Eq. (|30|) . this relation shows that H^ 1 (q, u>) vanishes at q = for any value of w: since it is an even 
function of q, this implies that it is proportional to q 2 in the small q limit for any lu. Thus, the condition (|A11|) 
resulting from the conservation law is satisfied. This result can be readily generalized to the case where a finite range 
particle-hole interaction is present. 



Appendix D: Expansion for small momenta 

A development up to second order in momentum q/kp is performed. 

1. Definitions 



We introduce the variables: 

k = k/kp 

X = q/kF 
£ = w/e kp 
7 = A/e kF 

X 2 2 

a = T - 7 

In the following, the reduced cutoff momentum will be denoted k c - We also define the functions: 



(Dl) 



e(K,7) = y (k 2 — l) 2 + 7 2 

g(n) = k 2 - 1 

P(k,7,0 = (e~.9) 2 x [{2ef-e 

R(K,y) = (e-. 9 ) 2 + 7 2 



It is useful to define the following integral: 



2 = 



dg 



Introducing the variables: 



the integral I becomes: 



f 2 

x = (e-g) - — + 7 
x\ = x(k c ) 

x 2 = x(k = 0) = (v/l+7 2 + l) 2 
f 2 

a = 7 c 

4 ' ^ 



dx 



i 2 2 

y + 7 



so that, defining the variable b = y/ 
if a > : Zox = 
if a < : Iqi = 



we have: 



■''2 



dx 



-L[ 2 dx J: L_ 

x 2 -b 2 2b J X1 \x-b x + b 



■'■■2 



dx 



c 2 + b 2 b 



arctan 



1 

21, 



In 



x — b 



x + b 



(D2) 
(D3) 
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Two approximations will be used to obtain explicit expressions for the integrals involved in the leading terms of 
the expansions: 

(1) n 2 dn ~ kcIk, corrcponding to an integrand peaked around the Fermi momentum (k = 1); 

(2) an integral can be neglected if its integrand is an odd function of k 2 — 1 (so that the integral would be exactly 
zero if its borns were ] — oo; +oo[ instead of [— 1; k c — 1]). 

We have checked the accuracy of these approximations. 

Note that some of the generalized Lindhard functions presented in this appendix contain a term which is a diverging 
integral, namely: [In |e — <?|]^!1q- However, it cancels with another term when they are injected in the RPA equation. 



Expanded expression of X^ 



X™ = - 



2. Expanded expressions of , X^ , X^_ h , 



and X° F 



4 



(2^) 2 e F 



e-g 



k=0 



4 + 3 ppp 



^X ppp + ( X 3 ) 



The integral T-pvv has the expression: 



(D4) 



4a 



+ 



pi 

x 



K =o X V 7 _ ~ 

«. i / 3 7 4 5 7 2 e 2 e 

K =o a 2 \ 4 16 32 

xC 2 -4 7 4 + 3 7 2 e 2 - V 

K=0 



P 2 



K=0 Ct 



■ *»v - it 



Expanded expression for X V J^ 



x pp = _ 



4 
(2n) 2 e F 



■ gW K K c =o + aI + x ppm + °(x 3 



where the integral I ppm has the expression: 



L ppm 



Ix i a 7 6 ,7 4 i 2 e 



10 



R 2 
1 

1 

P 

X 
IP 
1 



K=0 
K c 

K=0 



X - 



K=0 



re=0 



X - 



4^2 

8 7 6 

12 7 4 
I 2 " 



n 2 5^ 2 



16 32 



p2 

x 



K = 



p2 



16 7 6 - 16 7 4 £ 2 + 5 7 2 £ 4 - y 



k=o a 



- 3 7 2 f - \ 



(D5) 
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Expanded expression for X 1 ^ 1 



where the integral I p h m reads: 



X 



ph 



4 



7 2 I+^I phm + o( X 3 ) 



--phm 



= Ix 

1 

R 
1 



(27T)2 e p 

1 ^3 7 6 _7^ 7^ 
4 4 

6 7 2 



a 2 V £ 



P 2 
1 



P 3 
1 

P 

a; 
P 
1 



K=0 



«=0 



«=0 



12 7 4 
I 2 " 



8 7 6 

I 5 " 

12 7 4 

I 2 " 



1 



- - 5 7 ' 

3/ 5 7 4 7 2 g 2 

e + 4 8 



p2 

IE 



x (-16 7 6 + 12 7 4 £ 2 - 2 7 2 £ 4 ) 



p2 



1 /8 7 8 



18 7 6 + 8 7 4 C 2 - 7 2 C 4 



Expanded expression for -X^-jp 



X 



4 



G-F 



(2^) 2 e F 



x 7£ 



| + ^lGFm+o( X 3 ) 



where the integral Xgfui reads: 



- Ix ^( v ^2-T6" + 32 



+ 



K=0 



2f 

£ 4 

7 
i 2 



2 3 



1 



p2 

X 



K=0 



5^ _ 5f_ e 

£ 4 + 4£ 2 16 + 32 

4 



-4 7 4 + 3? 2 7 2 -| 



p2 



1 /2 7 6 9 7 4 „ 2, 2 

X - -75 f" + 2 7 2 ^ 2 - T 

=0 a V £ z 2 4 



(D6) 



(D7) 



3. Velocity of the Goldstone mode and pseudo-Goldstone mode 



In the case of a pure pairing interaction, the real part of the determinant of the RPA matrix (|26p diverges if: 



which reads 



1 — X^Upair — 



£ 2 2y 2 

^I+^I PPP + (X 3 ) = 



(D8) 



(D9) 
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In the limit of vanishing £, this gives: 



8 Zpp P (g->0) 
3 J(f->0) 



with 



so that: 



For 7 <C 1 we have: 



2(^0) = 

ippp(e^o) = 



■3 7 2 



i? 2 



27 4 



i? 3 



[1/i?] -> --g 
which leads to the Goldstone velocity: 



[l/i? 2 



Jl/g 2 ] 16, 4 [l/i? 3 ] 
7 [1/i?] 3 7 [1/i?] 



[1/i? 3 



r 4 

X 5 ^ 3 



q V3m 



(D10) 

(Dll) 
(D12) 

(D13) 

(D14) 

(D15) 



Let us now consider the pseudo-Goldstonc mode, which is present for the proton supcrfiuid when the Coulomb 
interaction is taken into account, but screened by the electron component. In the static limit for the electron free 
response function, we can express the Coulomb interaction between protons as: 

(D16) 



1 + U 0e v c 



To express the determinant of the RPA matrix, it is useful to introduce the reduced Coulomb interaction v c such 
that: 



47re 2 



q- x 

The corresponding reduced screened Coulomb interaction is: 

v c v a Vc 



v c Aire 

— ; v c = - 



k 2 



(1-2C 



o(x 3 )) 



(D17) 



(D18) 



X X + W 

where Q is related to the screening length of the electrons. The Coulomb interaction now enters the expression of the 
determinant D of the RPA matrix ([26]) as <=p-(l — + o(x 3 )), and we obtain: 



D 



€ I A I + X ( o-^PPP 1 g Q2 



^7 2 ^ppp) 



In the limit (x,0 —> (0,0), the determinant vanishes for 



8X D 



where the term 



31 



VcT 



9 1 

+ Q 2 



(D19) 



(D20) 



gives the Coulomb correction to the velocity of the Goldstone mode. 
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We are considering here the outer core of neutron stars. At even higher density, the P2 neutron pairing channel is expected 
to play a role. 



